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Abstract 

A brief introduction of the temporal filter based partially resolved numerical simulation/very large 
eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic 
subscale model and its advantages over the linear subscale eddy viscosity model are described. In 
addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three 
turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to 
illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 
single injector to further demonstrate the differences in the calculated flow fields resulting from the 
nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow 
generated in a single-element lean direct injection (LDI) combustor, the calculated result has 
demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, 
unsteady turbulent structures while using a relatively coarse grid. 

1.0 Introduction 

Many engineering applications of the computational fluid dynamics (CFD) for combustor flows need 
to capture the relatively large scales of unsteady, turbulent structures at both low and high Reynolds 
numbers to facilitate a higher fidelity analysis of the design. The conventional Reynolds-averaged Navier- 
Stokes (RANS) approach is known to be limited for this kind of task, because, by definition, the RANS 
solution can not contain the abovementioned flow information. Most recently, an approach called 
partially resolved numerical simulation (PRNS) (Refs. 1 and 2) has been developed for the simulation of 
very large eddy turbulence while using RANS type of grid resolution. Although PRNS is mainly aimed at 
the very large eddy simulation (VLES), hence notated as PRNS/VLES, it can easily be set to perform the 
large eddy simulation (LES) when the grid spacing reaches the traditional LES type of resolution. 
PRNS/VLES is based on the concept of temporal filtering to avoid the frequently overlooked issues that 
the traditional LES approach suffers from when a spatial filter is employed in conjunction with a spatially 
non-uniform grid. In PRNS/VLES, the larger time scales (or lower frequencies) of the turbulence are 
directly calculated, and the effects of the unresolved time scales of the turbulence are modeled by a 
dynamic equation system. The contents of both the resolved and unresolved turbulence are regulated by a 
“resolution control parameter” (RCP), which is related to the width of the temporal filter. 

The basic equations of PRNS/VLES and its subscale models are grid invariant, i.e., they do not have 
grid spacing as a parameter in their formal formulations. Therefore, it is possible to achieve a grid- 
independent numerical solution, and this is a major difference from the traditional LES approach and its 
variants. Another distinction is that PRNS/VLES enables us to perform unsteady RANS (URANS), 

VLES, LES, and even DNS in a unified way through the judicious selection of the value of RCP, in 
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conjunction with the employment of a grid spacing whose numerical resolution can consistently support 
the scale contents stipulated by the selected RCP. It should be noted that the PRNS/VLES approach is not 
a variant of the popular hybrid RANS/LES (Refs. 3, 4, 5, and 6). There is no enforced transition between 
the perceived RANS and LES domains. 

The subscale model is always one of the key issues in all turbulent numerical simulations. Some less 
than satisfactory large eddy simulations using the Smagorinsky eddy viscosity and its variants have been 
reported in the past (Refs. 7 and 8). Recently, we have also noticed that in a PRNS/VLES calculation of 
the turbulent pipe flow at a low Reynolds number Re x = 180 (based on the pipe radius and the skin 
friction velocity), the turbulent fluctuations were not sustainable over a long period of time and were 
eventually suppressed; and this can be attributed to the attempt of just using the eddy viscosity to account 
for all the subscale effects, even though the eddy viscosity is constructed from a k-z dynamic equation 
system. It is known that, in addition to the dissipative and diffusive effects accounted for through the eddy 
viscosity, the effects of anisotropy and rotation should also be included in the subscale model, especially 
when the simulation is a very large eddy simulation. To construct a more general relationship between the 
unresolved turbulent stresses and the resolved turbulent flow field, we have followed the analysis of the 
rational mechanics and obtained a general constitutive relationship (Refs. 9 and 10). This relationship 
indeed shows that, in addition to an eddy-viscosity term, there are several other terms representing the 
anisotropy and the rotation effects due to the interactions between the resolved and unresolved turbulence. 
They then introduce source terms in the momentum equation to sustain the turbulent fluctuations in the 
calculated flow field. The simulations presented in this paper will demonstrate this unique feature of the 
nonlinear subscale model. 

In the following, a brief description of the PRNS/VLES equations and the subscale model will first be 
presented, followed by an outline for the concurrent selections of RCP and the numerical grid. The results 
of simulations from three internal turbulent flows are then presented: the turbulent pipe flow, the turbulent 
swirling flow issued from a LM6000 single injector, and the flow generated in a single element lean direct 
injection (LDI) combustor. 

2.0 Temporal Filter Based PRNS/VLES Approach 

Using a homogeneous temporal filter G{t - 1 ' ) , the large time-scale turbulent variable c|) and its 
density-weighted variable <j) can be defined as 




( 1 ) 


where the integral is over the entire time domain and G satisfies the normalization condition: 

| G(t - t r )dt' = 1. There are several such temporal filters available, the simplest one is the top hat filter: 


G(t-t) = 


J 1/ A r , 

<N 

E-h 

< 

VI 

1 

1 o , 

otherwise 


( 2 ) 


where A T is the width of the top hat filter. Using this filter, the left part of Equation (1) becomes 

^ t+A T /2 

) = — f §(t ,Xj )dt’ . 

A t j 


( 3 ) 


t— Aj 12 
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Equation (3) reveals an unified feature of § and § , because they will become the exact Reynolds- 
averaged quantity and Favre-averaged quantity when A T —> go. On the other hand, they will become the 
instantaneous turbulent quantity as A T —> 0. For a finite they represent the quantities of large time- 
scale turbulence. 


2.1 Basic Equation 

Performing the filtering operation defined by Equation (1) on the Navier-Stokes equations, we obtain 
a set of exact, basic equations for the resolved, large time-scale turbulence ( § and <j>): 


f „ 2 „ ^ 


V 


(P^z ) i (p UfUj ^ , p i T'ijJ + 

(pc) t +(P S t S \i j P s kk -4i,i + 

P,f+(P«i),j = °> P = P R T, 

x . ij = p (UfUj - UiUj ) , cjj = p (Uje - u{e). 




2 \iSySy \lSfcfcSfj + Q , 

V 3 J 


(4) 

(5) 

( 6 ) 
(7) 


where sy = [uf j + Ujj )/2 . The symbols ( \ t and ( ) ti represent the temporal and spatial derivatives, 

respectively, p, u h T ' p , e, and g are the density, velocity, temperature, pressure, internal energy per unit 
mass, and the radiation rate, p and k are the viscosity and heat conductivity of the fluid. R is the gas 
constant. Xy and q t are the extra terms that were created by the process of temporally filtering the 

nonlinear Navier-Stokes equations. They represent the effects of unresolved, small time-scale turbulence. 
However, they are not in closed forms, and must be modeled. We call them the unresolved turbulent 
stresses and heat fluxes. 

It is clear that the basic Equations (4) to (7) which govern the large time-scale turbulence are not 
associated with the grid spacing, hence they are grid independent or grid invariant. This feature must be 
maintained in the development of the subscale models. It is important to note that in the temporally 
filtered turbulent field the filtering time scale A T is not the sole factor in determining the effective filtering 
length scale of the numerical simulation, and this gives leeway in arranging the grid spacing to carry out a 
physically meaningful numerical simulation. 


2.2 Nonlinear Subscale Model 

In order to obtain the solution for the large time-scale turbulence using PRNS/VLES Equations (4) to 
(6), we must model the unclosed terms defined in Equation (7): Xy = p {UfUj - UfUj ) , q t = p ( u^e - u{e ). 
There are several ways to model these terms. The more sophisticated method is to directly solve the 
transport equations for the unresolved turbulent stresses and heat fluxes, which can be derived from the 
Navier-Stokes equations (see Ref. 1). This method will require modeling of higher order unclosed terms 
appeared in the transport equation system of Xy and q t . A less complicated way is to use a general 
constitutive relationship between the unresolved turbulent stresses Xy and the strain rate of resolved large 
time-scale turbulence Sy , c by . This general constitutive relationship is then simplified according to the 
flow complexity by truncating the higher order nonlinear terms of Sy , (by . For example, the simplest 
form is just a linear relationship, which is the widely used subscale eddy viscosity model. Even at this 
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level, there are several ways to formulate the subscale eddy viscosity. The most popular one used in the 
traditional LES is the Smagorinsky model (Ref. 11) and its variations, which explicitly uses the local grid 
spacing A as the filtering length scale. A more sophisticated one is the one-equation model, such as the 
one proposed by Kim and Menon (Ref. 12), which solves the transport equation of the unresolved 
turbulent kinetic energy k , and using -Jk as the velocity scale, but still keeping the local grid spacing as 
the filtering length scale. 

In the context of PRNS/VLES, the subscale model is required to be a function of the width of the 
temporal filter A T , but independent of the local grid spacing. The constitutive relationship for the 
unresolved turbulent stresses is derived from the general constitutive relationship (Ref. 13) by invoking 
the realizability condition and the rapid distortion theory limit. The current subscale model contains 
linear, quadratic and cubic terms while the temporal filter width A T appears via the resolution control 
parameter RCP. 

2.2.1 Modeling of Unresolved Turbulent Stresses Ty 

The model proposed for PRNS/VLES is the following: 


T ij =“ 2 /l P~ (% ~^ijhk l^ + ^ij^kk 

— ^ 3/3 P ~~2^ik^kj ~ ®ik$kj ) 


( 8 ) 


-K 


_ £ 4 p 

2^5/5 P —sfk^*kj ^ik^km^mj ~ ^kl^lm^mk^ij + kl s (Sy ~ 5 fjSj c j c / 3) 


where, sy = ( u iyj + Ujj )/2 , <% = (u itj - ujj )/2 , II s = (s kk s mm - s M si k )/ 2 . The model coefficients 
Cju , A 3 and A 5 are constrained by the realizability condition and the rapid distortion theory limit. They 
are not arbitrary but formulated as (see Ref. 13): 




4.0 + As-U* 


M=- 


1-0 ~A}Cl 


k_ s * 

V £ 


0.5 + 1.5^J-n*5* 

E 2 


in which. 


k 2 

i-eqp— 

A,= 8 


_k 4 7S*S*+Q*Q* ’ 
P TT a 


( 9 ) 


4=V6coscp, cp = ^arccos(V61T*), W* - , 

U*=W) S*-=sy-hys kk 


(10) 

( 11 ) 


The coefficients fuf 3 , and f 5 are functions of Aj/T, i.e. the ratio of the temporal filter width to the global 
time scale of the turbulent flow of interest, where the global time scale T can be considered as the 
maximum integral time scale in the entire domain. These functions must have the following property: 
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fi 


( 12 ) 




This is because the unresolved turbulent stresses Xy must vanish when the filter width A T vanishes, and 
x y must approach the Reynolds-averaged stresses Ry as A T increases towards the global time scale T . In 
PRNS/VLES, A t /T< 1 , so we may take the following general expansion: 



q + q 




i = 1,3,5. 


(13) 


All Cq must be zero, because f) must be zero as A T goes to zero. If we retain the remaining two leading 
terms as an approximation and assume the first order derivative of f is zero at A T /T = 1 to reflect that 
fi reaches its maximum value of 1.0 at A T /T = 1 , then all /■ will have the same form: 




i = l, 3, 5 


(14) 


We call the ratio A r /T the resolution control parameter (RCP). It controls, at the governing equation level, 
the content of resolved large time-scale turbulence in the resolved field, which may, in theory, contain all 
length scales that have not been physically filtered. However, at the numerical solution level, the grid 
spacing will be imposed and its associated solution (say, the subscale eddy viscosity) together with the 
RCP will determine the effective content of the large-length scale turbulence. We will further discuss this 
parameter in Section 2.4.1. 

2.2.2 Nonlinear Interaction Between Resolved and Unresolved Turbulence 

It is important to identify the various physical interactions between the resolved and unresolved 
turbulent scales and to understand how these interactions are mimicked in the numerical simulation. In the 
momentum equation (4), these interactions are represented by the term Xyj , which is unclosed and must 
be modeled. In the traditional LES, this term is modeled via the isotropic eddy viscosity. Therefore, the 
effect of the interactions is accounted for only by a modification to the viscosity of the fluid. However, 
the real physical interactions are much more complex than this. In fact, the general constitutive 
relationship of Ty contains many more terms, in addition to a leading term that is related to an eddy 
viscosity. For example, the model of u y in Equation (8) has two parts: the linear term and the nonlinear 
(quadratic and cubic) terms. Each part plays different role in the momentum equation. The linear part 
leads to a term acting like an additional viscosity (called the subscale eddy viscosity); and the nonlinear 
part leads to terms acting like additional sources capable of promoting the resolved large scale turbulence. 
This can be clearly demonstrated by plugging the t y model into Equation (4), which yields 

( — 2 _ A 

+ 2(^+|a r )%--8y(^i + |i r )% +SJ , (15) 

V 3 fj 

where, 

vt = f\ q p— (16) 

8 
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( 17 ) 


!^3./3 P ^2 (^ik^kj ^ik ^kj ) | 




f _k 4 r„ „ „ 

~ j 2^5 /5 P"^3~ ^ik^frj ~ ^ik^kj "*" ^ik^km^mj ~ ^kl^lm^mk^ij ^ s ($ij ~ ^ij^kk /^) 


Apparently, the linear part of the model adds an additional subscale eddy viscosity p T (which is 
isotropic) to the viscosity of the fluid p ; and the nonlinear part provides a complex source term Sf , 
which accounts for the effects of anisotropy and rotation. 

We have noticed that, although different subscale eddy viscosity models have been used in different 
LES approaches, most of them have neglected or missed the source term Sf . Our study shows that this 
source term could become critically important for some flow simulations, especially for those at relatively 
low Reynolds numbers or flows with strong rotation, and this will be demonstrated in Section 3. 1.1. 2 and 
Section 3.2.2. 


2.2.3 Transport Equations for Subscale k - £ 

To complete the proposed model for x lh we need k and s, the unresolved turbulent kinetic energy and 
its dissipation rate. Their exact transport equations can be derived from the Navier-Stokes equations, and 
contain several higher order unclosed terms due to the temporal filtering operation. Here, we briefly 
describe the procedure of the derivation. The first step is to establish the transport equation for x ij9 
followed by a tracing operation to establish the equation for x u (which is 2 p k), and this leads to the 
transport equation for k. 

The exact transport equation for the unresolved turbulent stresses x z y = (pUiUj-pUfUj ) is 

x ij, t + (“* x ij ), k = D ij + +P ij~ P e ij ’ ( 1 8 ) 

where D (/ - , <J> ; y , P t] and ps y are the diffusion term, the pressure-strain correlation term, the 

production term, and the dissipation term, respectively. The following expressions indicate that 
all terms on the right hand side of the equation, except for the production term P ip are unclosed 
and must be modeled. 


A , 


-(p UiUjU k p ill'll j u k ) + \ 2[iu js ik -^d ik \iUjS i 


j ^\k$ik ~ ^ik\^S n 

V -3 


'-U,* 


■ (likU j + X jkUi ) k - ( puj 8 ik + pu, 8j k ~ puj 8 ik - pit , b jk J 


+ \2 \kUiS jk -—djkHUjSf, 


jk 2 ^jk^^n 


J Kk 


®ij=2ps ij -2ps ij 
kj — J; ~ tjk U i,k 


pe« = 


(sik u j,k + s jk u i,k ) _ \XS mm Sjj 2 \l(^S j k ltjj. + Sj k lij k j 


4. 

3 l 
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Now, if the diffusion term Dy is modeled by a gradient-type diffusion of Xy with the effective viscosity, 
p + \x T , then the trace of Equation (18) becomes 


d _ . d d 

— p/v H ow k = 

dt dxi dxf 


OXi 


Vi/ -pe 


(19) 


in which, 0 /7 has been neglected by ignoring the effect of compressibility on O y . 
The dissipation rate is defined by ps /z - / 2 , i.e., 



r - 2 - A 


r 2 ^ 

ps = 

2 p 5^- 5' z y — ~ P S mm $11 

v 3 ) 

— 

2p SijSij ~ — P S mm Su 
v d y 


( 20 ) 


A model transport equation for the dissipation rate 8 can be constructed by the analogy to Equation (19) 
as 


d _ d _ „ d 

PS H p Uf 8 = 




dXj 


dXi 


(fl + flr)— E 

OXi 


-C el T V S V j-C e2 £j- 


( 21 ) 


where C gl and C g2 are the model coefficients. We have adopted the commonly used values of C gl = 1.45 
and C g2 = 1.92 in the present work; while keeping in mind that they can be further constructed as 
functions of the local subscale turbulence quantities (Ref. 14). 

2.2.4 Modeling of Unresolved Turbulent Heat Fluxes q t 

A common practice in modeling the unresolved turbulent heat fluxes qi = p(w;e-W;S ) is to employ 
the following isotropic model: 

q ( = -k t ej (22) 


where k t is the eddy diffusivity for the heat, which is often modeled as k t =\i T / Pr r , where Pr r 
(about 0.9) is the turbulent Prandtl number. However, based on the analysis of the constitutive 
relationship (see Ref. 9 and 10), the simplest form that considers the effects of strain and rotation should 
be 


<?/ 


= -Kre i - K r ~(c\Sij + C2<% ) c 


(23) 


Where q and c 2 are yet to be determined coefficients. This more general model will result in 
modifications to both the diffusion term and the source term in Equation (5): 


/ __ 2 „ \ 

(pe) , + (P «/ e) . = ((k + k T )e,i) + I^kk + 2jl sySy - - j ds kk s u +Q + Sf 

’ 1 V 3 


(24) 


where the extra source term originated from the unresolved turbulent heat fluxes is 


Sf = 


K r-(ci%+c 2 <%)e 


(25) 
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So far, this source term has not been considered in most of the RANS and LES simulations. Furthermore, 
equations similar to Equations (23) to (25) can be applied to any filtered species equation. 

2.3 Time-Averaging Turbulent Quantities 

Since various filters can be used to define the resolved turbulent quantities, a practical yet somewhat 
overlooked question is how to conduct an apple-to-apple comparison between the calculated results and 
the experimental data of turbulent flows. 

Most of the experimental data are the “mean” values of turbulent variables (e.g., velocity, 
temperature, pressure, etc.); they are either the pure time-averaged values (for incompressible flows) or 
the density-weighted time-averaged values, i.e., the Favre-averaged values (for compressible flows): 

772 r/2 

( < t’) = ^- J H = } P^’ ( 26 ) 

-* 00 _ TI2 ->°° \ p ; _ T , 2 


where 4> represents an instantaneous turbulent quantity, it can be u t , u{Uj , or UfUjUj ^ , etc., T is the entire 

time domain. (c|>) is the pure time-averaged (Reynolds-averaging) quantity and [4>]is the density-weighted 

time-averaged (Favre-averaging) quantity. This type of experimental data is actually based on the 
assumption that the measured turbulent flows are statistically steady, or at least approximately steady. 
Otherwise, the experimental data of “mean” values must be redefined via the ensemble average of many 
repeated independent realizations of the (nominally) same experiment, i.e., 



where M is the total number of realizations and (]>M and (p^)^ are the values from each independent 

realization. The symbol { } represents the ensemble average. For statistically steady turbulent flows, the 
time average and the ensemble average become identical, i.e. 

« = <♦>. ^ = W ( 28 ) 

In summary, the experimental “mean” values for the compressible flows are either the density-weighted 
time average [())] or the density- weighted ensemble average {p4>} / {p} , expressed by Equations (26) and 
(27), respectively. 

Now, in PRNS/VFES, the temporally filtered values (|) , of turbulent variables are directly 
calculated. As a result, this enables an apple-to-apple comparison between the PRNS/VFES data and the 
experimental data. In Reference 1, relations between the PRNS/VFES data and the experimental data 
have been established for the first order, second order and higher order correlations. Furthermore, the 
relations for the Reynolds stresses and the scalar fluxes (which are second order moments formed by the 
zero-mean, fluctuating components of the velocities and scalars) can also be derived. Here we just present 
these relations under the statistically steady condition: 

h-] = (pM*)/(p)> H = (pe)/(p) (29) 

\ u i e \ = (p M / e ) / (p) > [“/“>] = (p“i“y)/(p)» [ u i u j u k] = (pUiUjMk)/{p), •••• (30) 
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Equation (29) indicates that the experimentally measured mean velocity and mean scalar (the terms 
on the left hand side of the equation) can be compared directly with the post-processed PRNS/VLES data 
(the terms on the right hand side). Even though Equation (30) cannot be used for direct comparison due to 
the unclosed terms on the right hand side, they are useful for establishing the following relationships for 
the Reynolds stresses and the scalar fluxes: 


(31) 

^=3/ +(«!>/( p). (32) 

where, 

R ij =(p( M * -[“f])(“y - [“y ]))/ (p) =[«/“y ] - [«/ ][“y ] » ( 33 > 

R i =^p(u i -[u i ])(e-[e]))/(p) = [u i e]-[u i ][e\, (34) 

Tij = (pu i u j )/(p)-(pu i )(pu j )/(p) 2 , (35) 

Ti = (pw, e)/(p) -(pM,-)(pe)/ (p) 2 . (36) 


The experimentally measured Reynolds stresses Ry and the scalar fluxes R t are determined according 
to Equations (33) and (34), respectively. Their PRNS/VLES counterparts, on the right hand side of 
Equations (31) and (32), consist of two components: the first component, Ty or 7), can be determined 
by using the resolved variables according to Equation (35) or (36); the second component, the 
unresolved turbulence stresses Xy or the scalar fluxes q t , will have to be obtained via the subscale models, 
Equations (8) and (23), discussed in the previous section. It is not difficult to recognize that the second 
component will evolve to the total Reynolds stresses and fluxes when the width of the temporal filter 
becomes sufficiently large, i.e., the contributions from both Ty and 7} should vanish. 


2.4 Guideline for Conducting the Simulation 

In the unfiltered turbulent flow field, the smallest time, length and velocity scales are the Kolmogorov 
micro-scales (Ref. 15). By invoking an analogy to the Kolmogorov micro-scales, the smallest time, length 
and velocity scales occurring in the filtered field are 


f 77 . 


V T + v 


L prns 


V 


J 


fl pr, 


(v T + v) 3 


Vprns =[(v r +v)c]4 


(37) 


Equation (37) indicates that, in addition to the viscosity of the fluid v, these smallest scales also depend 
on the mean values of subscale eddy viscosity v T and dissipation rate s , which are turbulent flow 
properties. Following Equation (37), the relationship between the length scale r\ pms and the time scale x pms 
is 


1 

r \prns ~ T 


(38) 
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Equation (38) indicates that for a fixed smallest time scale x prns (its minimum value should be A T in the 
filtered field), the smallest length scale x\ prns can automatically adjust itself to the change of the mean 
subscale eddy viscosity v T , or vice versa, i.e. the mean subscale eddy viscosity v T will adjust itself to 
the given local grid spacing which, in this scenario, is viewed as the local smallest length scale r \ prns . In 
this sense, there is some leeway in selecting the grid spacing for a given temporal filter width. Equation 
(38) also suggests that, for a given A r , when v T is no longer changing with respect to the reduction of 
grid spacing, the smallest length scale in the temporally filtered field will also not be affected by grid 
spacing reduction, i.e., a grid independent solution is reached. 

At the governing equation level, v T and s’ are functions of the resolution control parameter RCP 
(i.e., the temporal filter width). At the numerical solution level, their calculated values are functions of the 
RCP as well as the local grid spacing. The consistency conditions require that the resulted time scale x pms 
is of the same order of the temporal filter width A r , and the resulted length scale x\ prns can be of the same 
order of the grid spacing or larger. In practice, these two conditions can only be rigorously verified 
posterior, because the values of the time and length scales (' i prns , x\ prns ) will not be known until the 
simulation is completed. As a rule of thumb, a ‘successful’ simulation often is indicative of these two 
conditions being by and large satisfied. 

2.4.1 Resolution Control Parameter RCP 

The resolution control parameter RCP (i.e., the ratio A t !T) is used in PRNS/VLES to regulate the 
content of the time scales of the resolved turbulence. When RCP^> 1.0, all time scales of turbulence have 
been filtered (i.e., modeled), and the PRNS/VLES simulation becomes a RANS simulation, i.e., the 
directly calculated field is intrinsically the time mean, no turbulent fluctuation occurs in the directly 
calculated quantities. As the value of RCP decreases, the turbulent fluctuation becomes more pronounced 
in the directly calculated field, and the simulation moves towards VLES or LES. 

To carry out a very large eddy simulation, we need to choose a value of RCP from the outset. The 
consistency condition requires that x prns is of the same order of the temporal filter width A T . Therefore, by 
using Equation (37), we have 


Aj — 


l 


f V T + V^|2 



k 

8 


(39) 


This is based on the assumption that v T » v for a very large eddy simulation and the mean subscale eddy 
viscosity v T is of the order of k 2 /e , where k and s' are the mean value of subscale kinetic energy and 
its dissipation rate. As a result, the resolution control parameter RCP (A 7 / 7) is estimated according to 



k / s 

^ref I ^ ref 



(40) 


where T is the global time scale of interested turbulent flow and is expressed as a ratio of the reference 
turbulent kinetic energy to its dissipation rate. We have also assumed that s re f and s' are of the same 

order, since the dissipation mostly concentrates in the small scales (Ref. 15). Equation (40) suggests a 
way to guide the selection of the value of RCP. For example, RCP = 0.20 means that we intend to directly 
resolve those turbulence scales that are responsible for about 80 percent of the total turbulent kinetic 
energy while to model the rest unresolved turbulence scales that contain about 20 percent of the total 
turbulent kinetic energy. 

In the simulations presented in Section 3.0, the coefficientXAr/T) for the subscale model in Eq. ( 8 ) is 
set to be 0.3, 0.3333 and 0.4444 for the pipe flow, the LM6000 injector and the LDI combustor, 
respectively. According to Equation (14), the respective A t /T is about 0.16, 0.18 and 0.25, which means 
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that the PRNS/VLES simulation is to resolve the large scale turbulence that contains about 84 percent of 
the total turbulent kinetic energy for the pipe flow, 82 percent for the LM6000 injector and 75 percent for 
the LDI combustor. 

2.4.2 Arrangement of Numerical Grid Spacing and Grid Independent Solution 

Usually, a computational grid based on a priori knowledge of the flow field is generated subject to the 
constraint of the available computing resources. Typically, finer grid spacing is used in regions where the 
shear rates are higher and the turbulent fluctuations are stronger. Then we decide, as practically as 
possible, the extent to which the turbulence is to be resolved at the governing equation level and choose a 
value of RCP according to Equation (40). As discussed above, the consistency of the calculated results 
needs to be assessed posterior. In general, simulations with different levels of coarse grid spacing are all 
legitimate as long as the consistency conditions are satisfied. And of course, finer grid spacing is capable 
of revealing finer turbulent structures, but the statistical mean should be about the same. By refining the 
grid, it is possible to attain a grid independent solution; the criterion should be that for a given value of 
RCP the calculated mean values of x prns and r\ prns do not change when the grid spacing is further reduced. 

3.0 Numerical Results of PRNS/VLES Simulations 

The simulations of three internal flows will be presented, the first one is a turbulent pipe flow to 
illustrate the basic features of PRNS/VLES; the second one is a swirling turbulent flow associated with a 
GE LM6000 single injector to further demonstrate the advantage of the nonlinear model over the pure 
eddy viscosity model; and the third one is the flow in a single-element lean direct injection combustor. 

The last simulation was performed with a relatively coarse grid but for a quite complex geometry, and the 
dynamically important, unsteady turbulent flow structures have been revealed. 

3.1 Pipe Flows 

The turbulent pipe flow is one of the ideal benchmark flows for evaluating/validating numerical 
simulation approaches, because the geometry is simple and the flow is statistically homogeneous in the 
axial direction so that periodic boundary conditions can be applied to the inlet and the outlet to avoid the 
often complicated boundary-condition issues. In addition, some experimental data (Ref. 16) are available 
for comparison. 

Turbulent pipe flows at several different Reynolds numbers (from high to low) have been used to 
evaluate the fundamentals of the PRNS/VLES approach (Ref. 17). At each Reynolds number, the 
simulation was performed over a range of f (RCP) (from 0.2 to 1.0) to demonstrate the unified feature of 
the current approach. Also, two types of subscale models were used for flow simulation at each Reynolds 
number: one is a pure eddy viscosity model, and the other is a nonlinear model described by Equation (8). 
We have noticed that the nonlinear subscale model is critically important for the successful simulation 
of low Reynolds number pipe flow. Here we present some of the results at two Reynolds numbers: 

Re t = 180 and 3322 (based on the pipe radius and friction velocity). The value of f (RCP) was set to 0.3, 
which is a typical value for a very large eddy simulation. 

The computational domain has an aspect ratio of Length/Radius =10, where Radius = 0.06468 
m. Figure 1 shows the grid spacing in a cross section of the pipe. The grid is nearly orthogonal 
everywhere. The total number of the hexahedral elements is 906,750, more specifically, (156 x 29 + 39 x 
39) x 150. 
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Figure 1. — Grid spacing in a cross section of the pipe. 


This grid was used for flows at both Reynolds numbers 180 and 3322. The initial flow field was 
created by an arbitrary smooth profile plus random fluctuations. At the wall, a generalized wall function 
was imposed (Ref. 18). 

3.1.1 Reynolds Number Re r= 180 

Turbulent pipe flow at the low Reynolds number Re x = 180 (or Re = 7,000 based on the pipe 
diameter and centerline velocity) is relatively weak. In this section, we will present the numerical results, 
which demonstrate that for such a low Reynolds number turbulent flow, the pure eddy viscosity subscale 
model is unable to develop and sustain turbulent fluctuations in the numerical simulation. 

3. 1.1.1 Nonlinear Subscale Model 

The PRNS/VLES simulation using the nonlinear subscale model given by Equation (8) has 
successfully simulated the turbulent pipe flow. The transition from the initial disturbances to the fully 
developed, large scale turbulence can clearly be identified from the time histories of the velocity 
components (w, u and v). For a fully developed turbulent pipe flow, the turbulent fluctuations are stronger 
near the wall where the shear rates are higher. Near the centerline, the turbulent fluctuations are weaker. 
We have arranged 15 probes along a radius of the pipe to record the time history of the flow variables. 
Also, we have another 16 probes along the centerline to examine the statistical homogeneity in the axial 
direction. Here we only present the time histories of three velocity components at the probe 1, which is 
located at the center of the pipe. 

The time history of the velocity indicates that the initially induced random fluctuations are 
significantly damped during the first 600 time steps, then fluctuations quickly grow and evolve into a 
fully developed turbulence (See Figs. 2 and 3). 

The instantaneous contours of the axial velocity in a center plane (Fig. 4) and in a cross section (Fig. 
5) clearly demonstrate the motions of various turbulent scales. 
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Figure 3. — Time history of the u, v components. 
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Figure 4. — Contour of axial velocity in a center plane. 
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Figure 5. — Contour of axial velocity in a cross section. 


3.1. 1.2 Linear Subscale Model 

The PRNS/VLES simulation using the pure eddy viscosity subscale model is unable to develop the 
turbulent fluctuation into a fully developed turbulent pipe flow. The time histories of the three velocity 
components w, u and v at the center of the pipe indicate that the initially induced random fluctuations are 
quickly damped out and no new turbulent fluctuation emerges (see Figs. 6 and 7). 

Figures 8 and 9 are the contours of the axial velocity in a center plane and in a cross section of the 
pipe. These smooth contours do not show any turbulent fluctuating feature. 

To further study these phenomena, we have restarted a simulation using the eddy viscosity subscale 
model but starting from a fully developed turbulent flow field obtained in Section 3. 1.1.1. We have 
observed that this initial fully developed turbulence could not be sustained, in fact, it was quickly 
annihilated. Then we continued the simulation by turning on the nonlinear subscale model, and we found 
that the turbulent fluctuations quickly re-emerged and evolved into a fully developed turbulent pipe flow. 
The results of above simulations suggest that the pure eddy viscosity model is not suitable for 
PRNS/VLES simulation of relatively weak turbulent flow which has a low Reynolds number. 
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Figure 6. — Time history of the axial velocity w. 
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Figure 7. — Time history of the velocity components u, v. 
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Figure 8. — Contour of w component in a center plane. 


Figure 9. — Contour of w component in a cross section. 


3.1.2 Reynolds Number Re t = 3322 

Turbulent pipe flow at high Reynolds number Re t= 3322 (or Re = 150,000 based on the pipe 
diameter and centerline velocity) has quite strong turbulence. These strong turbulent fluctuations are easy 
to sustain in the numerical simulation, even with the use of pure eddy viscosity subscale model. We have 
performed simulations using/ (RCP) = 0.3 with both the nonlinear subscale model and the linear eddy 
viscosity model. In both cases, the initially induced random fluctuations were able to evolve into the fully 
developed turbulent fluctuations. Here we only present the numerical results from using the nonlinear 
subscale model. 

3. 1.2.1 Time History 

The time histories of the velocity components w, u, v and the subscale turbulent kinetic energy are 
recorded at 15 probes along a radius in a cross section at the mid-length of the pipe. Figures 10 to 15 
present the time histories of the axial velocity and the subscale turbulent kinetic energy k at the probes 1 , 
6, and 14, which correspond to the locations r/R = 0.0, 0.5, and 0.9743, or y + = 3322, 1661 and 85. 
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Probe 1 is at the center of the pipe, Probe 6 is at the midpoint between the center and the pipe wall, and 
Probe 14 is very close to the wall. The time histories at other 28 locations along the pipe radius and axis 
are also available. These time histories indicate that the initially induced disturbances of the velocity 
components are at first somewhat damped, followed by the emergence of new fluctuations, which are 
then quickly amplified and evolved into a fully developed turbulent pipe flow. The amplitude of the 
velocity fluctuations increases from the center towards the wall, so does the unresolved turbulent kinetic 
energy. 
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Figure 10 — Time history of w component at Probe 1 . 
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Figure 11 — Time history of w component at Probe 6. 
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Figure 12 — Time history of w component at Probe 14. 
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Figure 13. — Time history of k at Probe 1 . 
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Figure 14. — Time history of k at Probe 6. 
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Figure 15. — Time history of k at Probe 14. 
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Figure 1 6. — Contour of w component in a center plane. 
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Figure 17. — Contour of w component in a cross section. 


3. 1.2.2 Instantaneous Contour 

To examine the turbulent flow structures, we have plotted the instantaneous contours of various 
variables in a center plane and in a cross section. Here we present the instantaneous contours of the axial 
velocity component, the vorticity magnitude, the subscale turbulent kinetic energy and the effective 
subscale viscosity in Figures 16 to 23. The contours of the velocity component and the vorticity 
magnitude illustrate the features of the resolved large scale turbulent structures, whereas the subscale 
turbulent kinetic energy and the effective viscosity illustrate the highly non-uniform features of the 
unresolved scales. 
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Figure 18. — Contour of the vorticity in a center plane. 
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Figure 19. — Contour of the vorticity in a cross section. 
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Figure 20. — Contour of k in a center plane. 


Pipe flow, Re= 150,000, at time step 41,400 
PRNS, Non-linear model', Grid-B 
NCC-1.1.8, no smoothing, P-ref= 100,290 
2nd=0.0, 4th= 0.001, CFL=1.0, dt=2.5e-5 
Convergence order 3, or Max iter= 120 


Y 


mu 


n 


7.5000E-04 
7.0000E-04 
6.5000E-04 
6.0000E-04 
5.5000E-04 
5.0000E-04 
4.5000E-04 
4.0000E-04 
3.5000E-04 
3.0000E-04 
2.5000E-04 
2.0000E-04 
1 .5000E-04 
1 .0000E-04 
5.0000E-05 


if 


'S’® 



Pipe flow, Re= 150,000, at time step 41,400 
PRNS, NOn-linear model'. Grid-B 
NCC-1.1.8, no smoothing, P-ref= 100,290 
2nd=0.0, 4th=0.001, CFL=1.0, dt=2.5^5 
Convergence order 3, or Max iter^ 120 



turb_k 

n 4.5000E+00 
4.0000E+00 
h 3.5000E+00 
■ 3.0000E+00 
■ 2.5000E+00 
2.0000E+00 
i 1.5000E+00 
■ 1.0000E+00 
J 5.0000E-01 


Figure 21 . — Contour of k in a cross section. 


Pipe flow, Re= 150,000, at time step 41,400 
PRNS, NOn-linear model, Grid-B 
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Figure 22. — Contour of effective viscosity in a center plane. Figure 23. — Effective viscosity in a cross section. 
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3. 1.2.3 


Radial Profile 


The radial profiles of various flow variables (w, u, v, gauge pressure, vorticity magnitude, Mach 
number, subscale turbulent kinetic energy and dissipation rate, effective subscale eddy viscosity) are 
examined at different locations along the pipe axis z = 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.645 for the spatial 
development of the turbulent flow. Note that z is non-dimensionalized with the pipe radius. Figure 24 
shows the experimentally measured mean axial velocity (Ref. 16) together with the (calculated) 
instantaneous axial velocity at six downstream locations. Figure 25 presents the profiles diametrically. 
Figures 26 and 27 illustrate the radial profiles of the subscale turbulent kinetic energy and the effective 
viscosity, respectively. These results indicate that the simulated turbulent flow is fully developed and is 
statistically homogeneous along the pipe. 




Figure 24. — Radial profile of the axial velocity. Figure 25. — Radial profile of instantaneous axial velocity. 
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Figure 26. — Radial profile of instantaneous k. Figure 27. — Radial profile of instantaneous effective viscosity. 
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3. 1.2.4 Spectrum and Correlation Analysis 

Figure 28 and 29 present the power spectrum density (PSD) of the axial velocity component and its 
two-point (time) correlation at three locations (Probes 1, 6 and 14). The broadband feature of the PSD 
(i.e., more than two orders of energy variation from small scale to large scale) and the typical two-point 
correlation shapes (i.e., the correlation rapidly decreases as the time lag increases) indicate that the 
PRNS/VLES simulation with / (RCP) = 0.3 does mimic the statistical features of a fully developed 
turbulence. 


PSD of axial velocity w. 
Pipe flow, Re= 150,000 



HZ 


Figure 28. — Power spectrum density of the axial velocity 
at Probes 1 , 6 and 14. 



Figure 29. — Two-point time correlations of <ww> at Probes 1, 6 and 14. 
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Pipe flow, Re = 150,000 
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Figure 30. — Comparison of mean axial velocity between 
PRNS/VLES, URANS and experimental data. 


The experimental data of the mean axial velocity at Reynolds number of 145,700 is used to compare 
the time averaged axial velocity obtained from the current simulation at a Reynolds number of 150,000. 
In Figure 30, the result obtained from URANS using a standard k - 8 model and the wall function is 
included. The result obtained from PRNS/VLES is in reasonable agreement with the experimental data, 
while the URANS result exhibits significant under-prediction. 


3.2 Flow in a Single-Element LM6000 Injector 

LM6000 is a General Electric low NOx gas turbine. We have performed several types of numerical 
simulation of the flow issued from one of its fuel injector. A highly swirling jet is injected from a circular 
inlet. The inlet pressure is about six atmospheres and the inlet temperature is about 644 K. The inlet 
boundary condition is the specified mean profiles from the experiments. The combustor is a rectangular 
box. The Reynolds number based on the inlet axial velocity and the inlet jet diameter is about 3,200,000. 
Figure 31 depicts the computational domain and the numerical grid. The total number of the grid points is 
about 495,000. This is the only grid used in this study to perform all of the simulations including 
PRNS/VLES, URANS and RANS. Here we only present the results from the very large eddy simulation 
using the PRNS/VLES approach. 

The value of f (RCP) is set to 0.3333. A convective unsteady boundary condition (Ref. 19), similar to 
that proposed by Ferziger (Ref. 20) and Grinstein, et al (Ref. 21), is applied at the outlet boundary. The 
initial condition is the solution of a steady RANS simulation (see Figs. 32 and 33 for the axial velocity u 
and the vorticity magnitude in a center plane). Two subscale models (nonlinear and linear) have been 
used. The results indicate that the nonlinear model is very helpful in the simulation of swirling turbulent 
flow, which occurs commonly in the combustor system. 
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Figure 31. — Computational domain for LM6000 single 
injector flow simulation. 
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Figure 32. — Contour of the axial velocity (RANS). 
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Figure 33. — Contour of the vorticity magnitude (RANS). 


3.2.1 Nonlinear Subscale Model 

The time history of velocity components and gauge pressure are recorded at four locations along the 
centerline: x = 0.015, 0.05, 0.10 and 0.2 m. From which we may monitor the development of turbulent 
fluctuations. Here, the presented time histories of the velocity components are at x = 0.1 and 0.2 m, which 
are located before and after the rear stagnation point of the recirculation zone. Figures 34 and 35 indicate 
that the turbulent fluctuations are fully developed. 


NAS A/TM— 20 1 0-2 1 6323 


21 



Time histry at Probe 3 Time histry at Probe 4 





PRNS, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 70,000 time step 
starting from steady RANS solution 
Using convective BC at exit 




u 

1.05E+02 
9.56E+01 
8.60E+01 
7.64E+01 
6.67E+01 
5.71E+01 
4.75E+01 
3.78E+01 
2.82E+01 
1 .85E+01 
8.91E+00 


-7.26E-01 


-1.04E+01 


-2.00E+01 



Y 

■n 


PRNS, non-linear model 
2nd= 0.0, 4th=0.01 , cfl = 1.0, dt = 4.E-06 
0.5 M elements, at 70,000 time step 
starting from steady RANS solution 
Using convective BC at exit 


Vorticity Magnitude 

H 3.00E+04 

■ 2.77E+04 

■ 2.54E+04 
Xfl 2.31 E+04 

2.08E+04 
1.85E+04 
1.62E+04 
1.39E+04 
1.16E+04 
9.30E+03 
7.00E+03 
4.70E+03 
2.40E+03 
1.00E+02 



Figure 36. — Contour of instantaneous axial velocity. Figure 37. — Contour of instantaneous vorticity. 


The instantaneous contours of flow variables in a center plane are presented for the time step 60,000 
which is about 100 through flow times (TFT, defined as the ratio of the length of the combustor to the 
inlet centerline axial velocity). Figures 36 and 37 are the contours of the axial velocity and the vorticity 
magnitude, respectively. The recirculation zone, the massive separation and the shear layers are clearly 
visualized. 

3.2.2 Linear Subscale Model 

In the early stage of the simulation, the calculated turbulent structures look reasonable. However, they 
are not sustainable in the long run, and the calculation crashes after 18,000 time steps (about 36 times of 
the through flow time, TFT). 

Here, we present the instantaneous contours at two instances: the time step 10,000 (about 20 TFT) 
and the time step 18,000 (about 36 TFT). Around 20 TFT, the flow structures, shown in Figures 38, 40 
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and 42, look reasonable except for the subscale turbulent quantities. The subscale turbulent kinetic energy 
k (not shown) and the eddy viscosity \x T are too small (about two orders of magnitude smaller) when 
compared with its counterpart obtained from the nonlinear subscale model. Around 36 TFT, the flow 
structure near the outlet becomes unphysical (see Figs. 39, 41 and 43), e.g., large amount of high speed 
inflow occur at the exit. The subscale turbulent kinetic energy (not shown) and the eddy viscosity become 
even weaker and smaller. The simulation crashes soon after 36 TFT. These results suggest that the linear 
subscale k-s model does not work very well for highly swirling flows. 
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Figure 38. — Axial velocity at 10,000 time step. 
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Figure 39. — Axial velocity at 18,000 time step. 
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Figure 40. — Contour of vorticity at 10,000 time step. 


Figure 41 . — Contour of vorticity at 1 8,000 time step. 
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Figure 42. — Effective viscosity at 10,000 time step. 


Figure 43. — Effective viscosity at 18,000 time step. 


3.3 Flow in a Single-Element LDI Combustor 

The lean direct injection (LDI) injector is a liquid fuel injector developed to reduce aircraft emissions. 
Stable combustion is essentially completed within a short distance through rapid fuel and air mixing. This 
design also allows for many small fuel injectors integrated into modules facilitating different fuel staging 
strategies, such as the one shown in Figure 44. So far, experimental observations have not fully clarified 
the dynamics of the mixing and combustion processes occurring in these injectors, and numerical studies 
need to be conducted to achieve a better understanding of the underlying unsteady physics of the LDI 
combustor. 

A very large eddy simulation (PRNS/VLES) has been carried out for the non-reacting turbulent flow 
in a single-element LDI combustor as the first step towards the simulation of, for example, a 3x3 injector 
module. 

Figure 45 depicts the single-element LDI combustor geometry and its computational domain. The 
numerical grid is formed using hexahedral elements and the total number of elements is about 862,000, 
which is a relatively coarse grid used in a previous RANS simulation (Ref. 22). Embedded in this figure 
are the instantaneous iso-surface of zero axial velocity component colored by the subscale effective 
viscosity and six instantaneous stream lines originating from the inlet of the injector, then passing through 
the swirler and the convergent-divergent nozzle, finally entering the combustion chamber. 

In this study, we first carry out a steady RANS simulation to provide the initial flow field for the 
PRNS/VLES simulation. The inlet boundary condition is set by specifying the velocity, the density and 
the temperature based on the experimental data, the outlet boundary condition is an unsteady convective 
boundary condition (Ref. 19). The nonlinear subscale model is used since it has been proven to be very 
effective in our previous studies. The value of f (RCP) is set to 0.4444 since the available numerical grid 
spacing is very coarse. We have also carried out an unsteady RANS (URANS) simulation to compare 
with the PRNS/VLES simulation. Available experimental data is used to assess the simulation results. It 
is demonstrated that, even with a RANS type of grid spacing, the PRNS/VLES approach can successfully 
reveal the complex unsteady turbulent structures occurring in this single-element LDI combustor. These 
dynamically important flow structures include the precessing vortex core (PVC) and the vortex 
breakdown bubble (VBB). Good comparisons of velocity components with the experimental data are also 
demonstrated. 
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Figure 44. — Configuration of a 3 x 3 LDI injector module. 



Figure 45. — A snap shot of the flow and the grid spacing. 


3.3.1 Flow Structures 

In the following, the instantaneous contour snap shots and the iso-surfaces of some quantities are used 
to illustrate the flow structures. 

3.3. 1.1 Instantaneous Contour Plots 

Figure 46 is the contour plot of the axial velocity component in a center plane at time step 90,000. It 
shows that a strong recirculation zone is extended from the combustor dump plane deep into the upstream 
nozzle throat. Figure 47 is the contour plot of the subscale turbulent kinetic energy, which is high near the 
high shear regions. 


NAS A/TM— 20 1 0-2 1 6323 


25 



u 

B 9.14E+01 
8.17E+01 
7.20E+01 
6.24E+01 
5.27E+01 
4.30E+01 
3.33E+01 
^ 2.37E+01 
1.40E+01 

■ 4.34E+00 
-5.33E+00 
-1.50E+01 


PRNS (non-linear) of LDI 

2nd=0.0, 4th=0.01 , CFL=1 .0, dt=1 .Oe-6 

90,000 time steps 



turb_k 

S I .00E+02 

9.09E+01 
8.18E+01 
7.27E+01 
6.36E+01 
5.45E+01 
4.55E+01 
3.64E+01 
2.73E+01 
1 .82E+01 
9.09E+00 
0.00E+00 


PRNS (non-linear) of LDI 

2nd=0.0, 4th=0.01 , CFL=1 .0, dt=1.0e-6 

90,000 time steps 




Figure 46. — Contour of axial velocity in a center plane. 


Figure 47. — Contour of subscale kina center plane. 


H 


pg 

1.41E+04 

1.26E+04 

1.10E+04 

9.45E+03 

7.89E+03 

6.34E+03 

4.78E+03 

3.23E+03 

1.67E+03 

1.13E+02 

-1.44E+03 

-3.00E+03 


PRNS of LDI, dt=1.0e-6 
2nd=-0.0, 4th=0.01 , CFL=1.0 
90,000 time step 





Figure 48. — Side view of the PVC and the VBB. Figure 49. — Forty-five degree view of the PVC and the VBB. 

3.3.1.2 Instantaneous PVC and VBB 

The dominant flow structures in the LDI combustor can be best visualized via the iso-surface of the 
zero axial velocity and the iso-surface of a relatively low pressure. The iso-surface of the zero axial 
velocity is also known as the vortex breakdown bubble (VBB). The iso-surface of a sufficiently low 
pressure captures the precessing vortex core (PVC). Figures 48 and 49 are the snap shots taken from two 
different angles. Figure 48 is a side view and Figure 49 is a perspective view. In these figures, the dark 
blue region is a vortex core, which is formed near the nozzle throat and extends into the combustor 
chamber. This spiraling vortex rotates and breaks, it changes randomly in space and time. Embedded in 
these figures is an instantaneous stream line, which starts from the inlet of the injector and goes through a 
complex seemingly random path in the combustor chamber. This stream line spirals around the dark blue 
surface indicating that the dark blue region is indeed a vortex core. The light green surfaces are the iso- 
surfaces of the zero axial velocity. 

3.3.2 Profiles of Velocity Components 

The experimental data reported in Reference 23 is used to assess the current numerical results. The 
experimental data are mean values, but the numerical results are instantaneous values. 


NAS A/TM— 20 1 0-2 1 6323 


26 



3.3.2.1 


Axial Velocity Distribution Along the Centerline 

Figure 50 is a comparison of the axial velocity profile between the calculated instantaneous values 
and the experimental mean value. Figure 51 is an enlarged view near the dump plane of the combustor 
chamber. These figures show that the calculated values are in a reasonable agreement with the 
experimental data. It also supports the experimental observation that, near the dump plane, the turbulent 
fluctuations are quite large. 


3.3.2.2 Axial Velocity Distribution Along the y-Axis at Several Downstream Locations 

The distributions of the axial velocity along the y-axis in the cross section plane at several 
downstream locations x = 3-, 6-, 9-, 12-, 90- and 180-mm are presented for 10 different instants and 
compared with the experimental mean value. Figures 52, 53 and 54 clearly indicate that the turbulent 
fluctuations are quite large near the inlet of the combustion chamber and are quickly reduced towards 
downstream as shown in Figures 55, 56 and 57. In addition, the largest turbulent fluctuations are off the 
centerline, somewhere between the centerline and the wall. 
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Figure 52. — Axial velocity at x = 3 mm. 
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Figure 53. — Axial velocity at x = 6 mm. 
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Figure 56. — Axial velocity at x = 90 mm. 


Figure 57. — Axial velocity at x = 180 mm. 


3.3.2.3 Velocity Components w and v Along the y-Axis at Several Downstream Locations 

The distributions of the other two velocity components w and v along the y-axis in the cross section at 
several downstream locations x = 3-, 6-, 9-, 12-, 15- and 90-mm are shown for 10 different instants and 
compared with the experimental mean value in Figures 58 to 69. Again, these figures clearly indicate that 
the turbulent fluctuations are larger near the inlet of the combustor chamber and are quickly reduced 
downstream. The strongest turbulent fluctuations are found off the centerline, somewhere between the 
centerline and the wall. 
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Figure 58. — Velocity component w at x = 3 mm. 
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Figure 59. — Velocity component w at x = 6 mm. 



Figure 60. — Velocity component w at x = 9 mm. 



Figure 61 . — Velocity component w at x = 12 mm. 




Figure 62. — Velocity component w at x = 15 mm. Figure 63. — Velocity component w at x = 90 mm. 
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Figure 64. — Velocity component v at x = 3 mm. 
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Figure 65. — Velocity component v at x = 6 mm. 



Figure 66. — Velocity component v at x = 9 mm. Figure 67. — Velocity component v at x = 12 mm. 



Figure 68. — Velocity component v at x = 15 mm. Figure 69. — Velocity component v at x = 90 mm. 
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3.3.3 Comparisons Between RANS, URANS, PRNS/VLES and Experimental Data 

This section will compare some time averaged mean velocity profiles from PRNS/VLES simulation. 
The mean values are obtained by time averaging the PRNS/VLES simulation over the last 10,000 time 
steps. 


3.3.3. 1 Axial Mean Velocity Along the y-Axis at Several Downstream Locations 

The distributions of the axial mean velocity along the y-axis in the cross section plane at several 
downstream locations x = 3-, 6-, 9-, 12-, 90- and 180-mm are presented for RANS, URANS and 
PRNS/VLES and compared with the experimental mean values. Figures 70 to 75 clearly show that the 
time averaged axial velocity profile from PRNS/VLES simulation are much closer to the experimental 
data, especially in the region near the inlet of the combustor chamber where turbulent fluctuations are 
strong. 



Figure 70. — Velocity component U at x = 3 mm. Figure 71 . — Velocity component U at x = 6 mm. 



Figure 72. — Velocity component U at x = 9 mm. Figure 73. — Velocity component U at x = 12 mm. 
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Figure 74. — Velocity component U at x = 90 mm. 
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Figure 75. — Velocity component U at x = 180 mm. 



3.3.3.2 Mean Velocity Components w and v Along the y-Axis at Several Downstream Locations 

The distributions of the other two mean velocity components w and v along the y-axis in the cross 
section at several downstream locations x = 3-, 6-, 9-, 12-, 15- and 90-mm are compared with the 
experimental mean value in Figures 76 to 87. Again, these figures clearly indicate that the time averaged 
mean velocity profiles from PRNS/VLES simulation are much closer to the experimental data, especially 
in the region near the inlet of the combustor chamber where turbulent fluctuations are strong. 
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Figure 78. — Velocity component V at x = 9 mm. 


Figure 79. — Velocity component V at x = 12 mm. 
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Figure 82. — Velocity component w at x = 3 mm. 
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Figure 83. — Velocity component w at x = 6 mm. 
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4.0 Conclusions 

The basic equations of the PRNS/VLES approach for the large or very large simulation are presented. 
They are based on the temporal filtering with a constant filter width. Consequently, they are grid-spacing 
independent or grid invariant. This feature allows the possibility of achieving a grid independent solution. 

The nonlinear subscale model is mathematically and physically more sound than the pure eddy 
viscosity model. The advantage of the nonlinear subscale model over the eddy viscosity model has been 
demonstrated in the simulations of the turbulent pipe flow at low Reynolds number Re x = 1 80 and the 
highly swirling flow issued from a LM6000 single injector. In both cases, the pure eddy viscosity model 
does not lead to a sustainable and physically meaningful solution. 

The simulations of the single-element LDI injector flow using the nonlinear subscale model has 
demonstrated that the PRNS/VLES approach can capture the dynamically important unsteady turbulent 
structures even with a grid spacing typically used in the RANS calculation. This is particularly 
encouraging, because the capability of revealing unsteady turbulent flow structures with a coarse grid is 
very much desired for the practical engineering application. 
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